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Abstract. - Multi- fractal model for hydrodynamic fully-developed turbulence (FDT) has been 
used to provide a detailed structure for the critical exponent a describing the scaling form of energy 
(or enstrophy) dissipation rate e (or r) that appears to exhibit an interesting universality covering 
radically different hydrodynamic FDT systems. This result also appears to provide a consistent 
framework for classification of dissipation field into critical, subcritical and supercritical cases. Some 
FDT problems that exemplify these cases are discussed. 

PACS numbers: 47.27.Ak, 47.27.Gs 

Introduction. - Small-scale structure in three-dimensional (3D) incompressible fully-developed 
turbulence (FDT), following Kolmogorov's [1] epoch-making work, is believed to possess, in the large 
Reynolds number (R ^ oo) limit, a certain universality in its scaling properties. This universal 
scaling behavior depends only on symmetries and conservation laws of the system and is unaffected 
by the large-scale flow structure. 

This is reminiscent of the scaling behavior near the critical point where many diverse systems 
show a striking similarity in their behavior [2,3]. Critical phenomena had a theoretical breakthrough 
in the renormalization group (RG) [4,5] which was the culmination of the ideas of scaling and 
universality. The formal application of RG procedure was attempted for the FDT problem [6-9]. 
The goal is to determine the critical exponents (like those associated with correlation length in 
critical lattice spin systems) that are intrinsic features of the FDT dynamics and not artifacts 
of the large-scale stirring mechanisms and hence unify radically different FDT systems near their 
critical points. 

Spatial intermittency is a common feature of FDT and implies that turbulence activity at small 
scales is not distributed uniformly throughout space. This leads to a violation of an assumption in 
the Kolmogorov [1] theory that the statistical quantities show no dependence in the inertial range 
L ^ £ ^ r] on the large scale L (where the particular external stirring mechanisms generating FDT 
become influential) and the Kolmogorov microscale r] = [y^ je)^!'^ (where the viscous effects become 
important). Thus, if one views the Kolmogorov [1] theory as a mean field theory [10], the spatial 
intermittency aspects will be expected to define at least one additional universal scaling exponent 
[11]. Spatial intermittency effects can be very conveniently imagined to be related to the fractal 
aspects of the geometry of FDT [12]. The mean energy dissipation field e may then be assumed, in 
a first approximation, to be a homogeneous fractal [13], and more generally, a multi-fractal [14-17]. 
The latter idea has received experimental support [18]. 

In the i? =^ oo limit (which corresponds to the critical point for FDT) infinitely many length 
scales become important. So, by analogy with critical phenomena, one may expect many macro- 
scopic details to become irrelevant and the critical behaviors of radically different FDT systems to 
exhibit some universality. Further, this universality may be expected to be strongly rooted in the 
self-similarity of the system at the critical point just as that for a finite critical ferromagnet system 
via the multi- fractal structure of the e-field [19]. 
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One of the fundamental dynamical assumptions underlying FDT theory is that e remains finite 
in the i? =^ oo limit (see [20] for experimental support for this assumption [21] - the infinite R 
limit here is taken to imply =^ 0, see (5) below). This implies that if one identifies e as the order 
parameter, then the critical exponent a defined by 



(1) 



is 



(7 = 0. 



(2) 



This aspect has also been validated by the multi- fractal model for the e-field [23]. The purpose 
of this letter is to report that the multi-fractal model indeed provides a detailed structure for the 
critical exponent a that appears to exhibit an interesting universality covering radically different 
FDT systems. 

The Multi-fractal Model for Classical 3D FDT. - In the multi-fractal model for the 
classical 3D FDT system, one stipulates that the fine-scale regime of FDT possesses a range of 
scaling exponents cc G / = [ctmin, "max]- Each a e I has the support set S{a) C of fractal 
dimension /(a) [24] such that, as ^ ^ 0, the velocity increment has the scaling behavior 



5vie) ~ r. 



(3) 



The sets S{a) are nested so that S{a') C S{a), for a' < a. 

On extrapolating this multi-fractal scaling in the inertial range down to the Kolmogorov mi- 
croscale rj [25] , the latter is found to exhibit the scaling behavior [26] 



(I) ~ (i?i) ^ 



where -Ri is a mean Reynolds number. 



Ri 



(eL4)V3 



(4) 



(5) 



e being the mean energy dissipation rate. 

The moments of the velocity-gradient distribution are then given by 



dv 



dx 



j {Rx)-^^'^ 



j^[pa-p+3-/(a)] 



diJ,{a) 



(6) 



In order to determine the scaling behavior of A^^ , first note that the dominant contribution to the 
integral in (6), in the limit oo, corresponds via method of steepest descent to 



_d_ 
da 



pa-p + 3- f{a) 
1 + a 



0. 



(7) 



Next, the sums of the moments of the total energy dissipation E{i) e(£)i'^ occuring in N{i) 
boxes of size i covering the support of the measure e exhibit the following asymptotic scaling 
behavior [27] 

i=l 
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Dq being the generalized fractal dimension (GFD) of the e- field [28]. The dominant contribution 
to the integral in (8), in the limit £ ^ 0, corresponds again via method of steepest descent to 

(3a* + 2)q - /(a,) = {q - l)Dq (9a) 

where, 

The coincidence of the values of a* given by (7) and (9b) for which the integrands in (6) and 
(8) become extremum is insured by assuming the Kolmogorov refined similarity hypothesis [18] in 
the Kolmogorov microscale regime. 

Combining (8) and (9), we obtain - 

Dq(p~S)-5p+9 

A!^^ ~ (i?i) ^0+1 , Ri^oo (10) 



where Q is the root of 

o = £ 

Dq + 1 



Q= \^'r' - (11) 



Using (11), (10) may be rewritten as 

~ (i?i)^^'\ Ri^oo (12a) 

where, 

7W^-(p-3Q). (12b) 
Thus, the mean energy dissipation i^^2^^ has the following scaling behavior - 

i^a'^^ ~ {RiT\ Ri^oo (13a) 

where, 

ai = 3(Q - 1). (13b) 
On the other hand, corresponding to = 2, (11) yields 

Q = 1 (14) 
in agreement with [23]. Using (14), (13b) leads to 

ai = (15) 

and (13a) then leads to 

uA^2^ ~ {Rif ~ constant, ^ oo (16) 

validating the inviscid dissipation of energy in classical 3D FDT! Further, (14) implies that the 
e-field in classical 3D FDT has the GFD Dq equal to the information entropy dimension Di . 

Applications to Other FDT Systems. - Let us consider next the classical two-dimensional 
(2D) incompressible FDT system [29]. Following the above procedure, the scaling behavior of the 
mean enstrophy dissipation uA':^^ (which is the relevant physical quantity in the enstrophy cascade 
of 2D FDT) can be calculated from the results in [30] - 

~(i?2)^^ i?2^oo (17a) 
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where, 

= 3(g - 1), (17b) 



and 

(fL6)V3 



R2 ~ ^ — (18) 



T being the mean enstrophy dissipation rate, and Q is the root of - 

^ Dq + a ■ ^'^^ 

(72 is identical to the exponent ai (13b) for the mean energy dissipation in 3D FDT! Corresponding 
to p = 2, (19) yields 

Q = l. (20) 

Using (20), (17b) leads to 

a2 = (21) 



uA^^^ ~ {R2f ~ constant, i?2 oo (22) 



and (17a) then leads to 

establishing the inviscid dissipation of enstrophy in the enstrophy cascade of classical 2D FDT. 

Let us consider next 3D compressible (barotropic fluid) FDT [31]. The scaling behavior of the 
mean kinetic energy dissipation ij,A2^^ can be calculated from the results in [31] - 

hA'^^ ~ {Ri^^, Ri^oo (23) 

where, /37-l\ 

as^[-^){Q-l) (24) 



and Q is the root of 

7 being the ratio of specific heats of the fluid. (24) exhibits the universality but with a different 
amplitude. The latter reflects the residual effect of variations in the cascade physics involving the 
acoustic wave dynamics in the compressible case [32]. Corresponding to p = 2, (25) yields 



Using (26), (24) yields 
and (23) leads to 



Q = 1. (26) 
CT3 = (27) 



hA''2^ ~ (i?i)" ~ constant, i?i ^ oo (28) 



establishing the inviscid dissipation of kinetic energy in 3D compressible FDT. 

These results imply that energy cascades of 3D and the enstrophy cascade of 2D FDT are indeed 
examples of the critical dissipation field which - 

• has the GFD equal to the information entropy dimension Di , 

• has a finite mean value in the inviscid limit. 
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This suggests further that we may define a suhcritical/ supercritical dissipation field as one which - 



• has the GFD less/ greater thsM the information entropy dimension 

• has an infinite/vanishing mean value in the inviscid limit; 

As an example of a suhcritical dissipation field, consider the 2D quasi-geostrophic FDT system 
in which the baroclinic effects are produced by the deformed free-surface of the ocean [30] . The 
scaling behavior of the mean enstrophy dissipation vA)^ can again be calculated from the results 
in [30] - 



vA. 



(2) 



{R2 



R2 



OO 



where, 



and Q is the root of 



(73 = 3(Q - 1) 



Q = 



3Dq + lip - 6 



(29a) 



(29b) 



(30) 



as is again identical to the corresponding exponents for the classical 3D and 2D FDT! Corresponding 
to p = 2, (30) yields, 



Q = l + 



Using (31), (29b) leads to 



and (29a) then leads to 



3Dq + 12 
(73 >0 



> 1. 



vAu^^ ^ OO, i?2 ^ OO- 



(31) 

(32) 
(33) 



So, the enstrophy dissipation in the enstrophy cascade for the 2D quasi-geostrophic FDT system 
diverges [36] in the inviscid limit in accordance with the first of the above stipulations for the 
suhcritical dissipation field. Further, (31) implies that the enstrophy dissipation field has the GFD 
Dq less than the information entropy dimension Di {Dq can be shown [37] to be monotonically 
decreasing with Q), in accordance with the second of the above stipulations for the suhcritical 
dissipation field. 

The Probability Density Functions of the Flow-Varaible Gradients. - The universality 
of the spatial intermittency aspects in FDT systems, on the other hand, would imply that the 
criticality, suhcriticality or supercriticality of the dissipation field should be traceable to properties 
like the probability density function (PDF) in the zero intermittency limit of the FDT system in 
question. 

The PDF of the velocity gradient in the zero intermittency limit for the classical 3D FDT is 
given by (Frisch and She [38]) 



ns) - ( ^ 



1/3 



y2/3| 



,|4/3' 



(34) 



where. 



s ~ 



dv 
dx 



and vq is the velocity increment characterizing large scales. 
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The PDF of the vorticity gradient in the zero intermittency Umit for the classical 2D FDT is 
given by (Shivamoggi [30]) 

"j^2/3|^|4/3' 



Pir) - ( rr 



1/3 



(35) 



where, 



dx- 



The PDF of the velocity gradient in the zero intermittency limit for 4D compressible FDT is 
given by (Shivamoggi [31]) 



Pis) 



1/3 



(36) 



where is a reference kinematic viscosity. 

The identity of (34)- (36) appears to confirm that energy cascades of 3D and the enstrophy 
cascade of 2D FDT are examples of the critical dissipation field. 

On the other hand, the PDF of the vorticity gradient in the zero intermittency limit for the 2D 
geostrophic FDT is given by (Shivamoggi [30]) 



5/11 



j,10/ll|^|12/ll 



(37) 



which is more non-Gaussian than (35) and is consistent with the fact that 2D geostrophic FDT is 
an example of a subcritical dissipation field. 

Discussion. - According to the multi-fractal model, the order parameter e (or r) in a variety 
of hydrodynamic FDT cases therefore has a scaling form 



with the apparently universal critical exponent a given by 

a = a{Q- 1) 



(38a) 



(38b) 



Dq being the GFD of the dissipation field. Q ^ I corresponds to supercritical, critical (of which clas- 
sical 3D and enstrophy cascade of 2D FDT are examples) , and subcritical (of which 2D geostrophic 
FDT is an example) dissipation-field cases. On first thought, the universality of critical exponent 
(T would appear to be rather far-fetched because the various cases of FDT are far too diverse to 
be characterized within a universality class. But, as speculated in [19], by focusing on the critical 
point {R =^ oo) behavior, it is apparently possible to do just that. 

I am thankful to Professors John Cannon, Greg Eyink, Grisha Falkovich, Siegfried Grossmann, 
Mike Johnson, David Kaup, Evgeny Novikov, Katepalli Sreenivasan, Henri Tasso, Angelo Vulpiani 
and Victor Yakhot for their helpful comments and advice. 
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